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Abstract 

In  this  paper,  we  are  interested  in  the  multiobjective  evaluation  of  the  schedule  performance  in  the  flexible  job  shops.  The  Flexible  Job 
Shop  Scheduling  Problem  (FJSP)  is  known  in  the  litteratue  as  one  of  the  hardest  combinatorial  optimization  problems  and  presents  many 
objectives  to  be  optimized.  In  this  way,  we  aim  to  determine  a  set  of  lower  bounds  for  certain  criteria  which  will  be  able  to  characterize 
the  feasible  solutions  of  such  a  problem.  The  studied  criteria  are  the  following:  the  makespan,  the  workload  of  the  critical  machine,  and 
the  total  workload  of  all  the  machines.  Our  study  relates  to  the  determination  of  a  practical  method  in  order  to  evaluate  the  representative 
performance  of  the  production  system. 
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T.  Introduction 

THE  flexible  job  shop  scheduling  problem  is  a  problem  of  planning  and  organization  of  a  set  of  tasks  to  be  performed 
on  a  set  of  resources  with  variable  performances  [1],  [2],  [3],  [4],  [5],  Tn  the  literature,  the  authors  generally  consider 
two  steps  in  its  resolution  [4],  The  first  one  is  to  assign  the  various  tasks  to  the  suitable  resources.  The  second  step  is 
the  sequencing  of  the  tasks  and  the  computation  of  the  starting  times  by  taking  account  of  the  various  constraints  of 
precedence  and  resources.  Nevertheless,  Dauzere-Peres  et  al  have  recently  proposed  an  interesting  method  to  solve  this 
problem  by  considering  the  two  steps  at  the  same  time  [5] .  Morover,  in  [6] ,  Dauzere-Peres  et  al  have  given  an  extension 
of  this  problem  in  which  tasks  can  be  performed  by  several  resources  at  the  same  time.  Many  criteria  can  be  considered 
in  the  resolution  of  such  a  problem.  Mainly,  two  objectives  could  be  distinguished.  The  first  one  is  to  balance  the 
workloads  of  the  machines,  the  second  one  is  to  organize  the  various  tasks  in  minimizing  the  overall  completion  times. 
This  variety  of  criteria  induces  additionnal  difficulties  related  to  the  evaluation  of  the  feasible  solutions  as  well  as  the 
comparison  between  the  resolution  methods. 

Tn  this  way,  we  aim  in  this  paper  to  contribute  in  the  reduction  of  such  difficulties  by  proposing  a  set  of  lower  bounds 
for  some  representative  criteria  of  such  a  problem.  Tn  addition,  we  propose  to  undervalue  the  criteria  related  to  the 
workload  of  the  resources.  The  second  part  of  this  article  presents  the  specificities  of  the  flexible  job  shops  and  gives  the 
mathematical  formulation  used  to  deal  with  such  a  problem.  Thus  in  section  TTT,  we  describe  the  various  steps  followed 
in  the  calculation  of  the  lower  bounds  proposed.  The  last  part  will  be  devoted  to  the  presentation  of  some  results  and 
some  conclusions  concerning  this  research  work. 

TT.  Problem  Formulation 

The  problem  is  to  organize  the  execution  of  N  jobs  on  M  machines.  The  set  of  machines  is  noted  U.  Each  job 
Jj  represents  a  number  of  rij  non  preemptable  ordered  operations  (precedence  constraint).  The  execution  of  the  ith 
operation  of  job  Jj  (noted  (),:l  )  requires  one  resource  or  machine  selected  from  a  set  of  available  machines.  The 
assignment  of  the  operation  Ojj  to  the  machine  M],  entails  the  occupation  of  this  machine  during  a  processing  time 
called  dijk-  Thus,  to  each  FJSP  (Flexible  Job  Shop  Scheduling  Problem),  we  can  associate  a  table  D  of  processing 
times  such  that:  D  =  {d,  j  /,.  £  IN*  |  1  <  j  <  TV;  1  <  i  <  rij ;  1:  <  k  <  M}. 

Tn  this  problem,  we  make  the  following  hypotheses: 

-  all  machines  are  available  at  t  =  0  and  each  job  Jj  can  be  started  at  t  =  fj , 

-  at  a  given  time,  a  machine  can  only  execute  one  operation:  it  becomes  available  to  other  operations  once  the 
operation  which  is  currently  assigned  to  is  completed  (resource  constraints), 

-  to  each  operation  Ojj,  we  associate  an  earliest  starting  time  r, j  calculated  by  the  following  formula: 
r\  j  =  r j  VI  <  j  <  N  ,  and  rj+ij  =  +  li  j  V  1  <  i  <  n j  —  1 ,  VI  <j<  N. 

where  ji  j  =  min*  VI  <  i  <  rij  —  1  ,V  1  <j<  N. 

The  FJSP.s  present  two  difficulties.  The  first  one  is  to  assign  each  operation  Ojj  to  a  machine  M*  (selected  from  the 

set  U).  The  second  one  is  the  computation  of  the  starting  time  t, y  and  the  completion  time  t fjj  of  each  operation  Oj  j  . 

The  considered  objective  is  to  minimize  the  following  criteria: 

-  the  makespan: 
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Cn  =  max  {//„  .  } 


(1) 


-  the  workload  of  the  most  loaded  machine: 


Cr-2  =  max  {114} 

k 


(2) 


where  114  is  the  workload  of  M]~, 

-  the  total  workload  of  the  machines: 


Cr3  =  V  (3) 

k 

Def  inition  1 :  To  each  repartition  of  the  operations  on  the  resources  set,  we  associate  an  assignment.  Each  assignment 
is  characterized  by  a  set  S  such  that:  S  =  ;c  <E  {0, 1 }  j  1  <j<  N,  1  <  i  <  nj  ,1  <  k  <  M } .  S'; y  *  =  1  if  Of  j  is  assigned 
to  M*  else  Si  jj.  =  0. 

Def  inition  2:  E  is  a  numerical  function  defined  as  follows:  E  (x)  =  x  if  x  is  integer  else  E  (x)  =  E  (x)  +  1,  where 
E  (x)  is  the  integer  part  of  x. 


Til.  New  Lower  Bounds 

The  problem  of  the  lower  bounds  has  been  considered  in  the  literature  for  many  scheduling  problems,  in  particular, 
the  one- machine  problem  [7],  the  parallel  machines  problrm  [S],  [9],  the  hybrid  flow-shop  problem  [11]  and  the  job- 
shop  problem  [12].  Generally,  the  methods  suggested  are  based  on  the  constraint  relaxation  (preemption  of  the  tasks, 
disjunctive  constraint  on  the  resources...)  in  order  to  estimate  the  makespan  of  the  optimal  schedule.  In  this  paper, 
we  generalize  some  results  proposed  in  the  literature  for  the  parallel  machines  problem  [9]  and  we  propose  others  new 
considerations  based  on  the  evaluation  of  the  cost  of  assigning  a  certain  number  of  tasks  to  certain  resources. 

A.  Lower  Bound  for  the  Workload  of  the  Machines 

The  total  workload  is  equal  to  the  sum  of  all  the  processing  times  of  all  the  operations  carried  out  by  the  set  of  the 
machines  according  to  the  chosen  assignment.  In  addition,  for  each  operation  0;,j,  the  processing  time  is  superior  or 
equal  to  y,  •,  therefore: 


C?’3  -  US T'iJ 

j  i 

B.  Lower  Bound  for  the  Workload  of  the  Critical  Machine 
Lemma  1:  E  I  '  J  ,|,)  1  is  a  lower  bound  of  Cr-2. 

Proof:  Trivial,  the  workload  of  the  critical  machine  is  higher  than  the  mean  of  the  workloads.  ■ 

Now,  let  consider  the  operations  carried  out  by  the  set  of  the  machines.  N  is  the  average  of  the  operation  numbers 
carried  out  by  one  machine  (N  =  -g.  £4  nj).  Then,  we  have  at  least  one  machine  Mj~0  that  will  perform  at  least  A40 

operations  such  that:  A40  >  N  =  E  (^N ^ .  Thus,  to  each  machine  M*,  we  associate  the  N  shortest  operations  that  it 
can  perform.  D}„  y  is  the  sum  of  the  processing  times  of  these  associated  operations  and  Sj  y  represents  the  minimal 
value  of  these  sums  when  varying  k: 

S,  a  =  min  (  D,  ) 

fto.A  l<k<M\  k’A  J 

The  critical  machine  will  have  a  workload  more  important  than  n  .  then,  we  obtain  the  following  lemma. 

Lemma  2:  Sko  y  is  a  lower  bound  of  CV2 . 
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Definition  3:  a  (A\  Nr)  is  a  binary  variable  which  is  true  if  the  machine  M carries  out  Nr  operations  among  the  NT 
operations  (NT  is  the  total  number  of  the  operations  and  Nr  <  NT). 

Let  us  suppose  that  a  (k,  Nf)  =  1,  then  the  machine  M} .  will  work  during  a  time  at  least  equal  to  D*  y;  (Dk.Ni  is 
the  sum  of  the  Nf  shortest  processing  times  of  the  operations  that  we  can  perform  on  Mk).  Tn  this  case,  the  remainder 
of  the  operations  (the  NT  —  Nr  operations)  must  be  carried  out  on  Uk  =  U  —  {JUT*}.  This  obligation  can  be  regarded 
as  a  scheduling  problem  of  (NT  —  Nf)  operations  on  the  (M  —  1)  machines  of  Uk .  Therefore,  for  such  a  problem,  we 
can  apply  the  result  of  lemma  1  to  undervalue  the  workload  of  the  most  loaded  machine  of  lJk .  Thus,  this  workload 
is  superior  or  equal  to  ak,Ni  =  E  ^  w^K:r<:  ^k,Ni  is  the  sum  of  the  (NT  —  Nr)  smallest  values  of  7)7  such  that 

7)  •  =  minj<A<M  (d%.j,h)-  Then,  we  obtain  the  following  lemma. 
h^k 

Lemma  3:  If  a  (k.  Nr)  =  1 ,  then  there  exists  a  machine  for  which  the  workload  is  at  least  equal  to  /3k  Nl  such  that: 

Pk,Ni  =  max  (Dk.Nh  £**, Nl) 

Lemma  j.  iL:(l  v  —  mini<A<M  ( min^ ,  (Jk.Ni  »  is  a  lower  bound  of  Cr 2  and  it  is  superior  to  6ko  N . 

Proof:  Obvious,  according  to  lemma  3  and  the  definition  of  A7,  ■ 


Using  lemmas  1  and  4,  we  deduce  the  following  lemma. 
Lemma  5: 


,  .  (  ::  ( xZj  xZi  ii.j 

Cr-2  >  max  - )  Nko,i v 


C.  Lower  Bound,  for  the  Make-span 
Lemma  6: 

Cr  1  >  max 

j 

Proof:  Obvious,  because  of  precedence  constraints  on  the  operations  over  the  jobs. 
Lemma,  1:  Let  RM  the  sum  of  the  M  smallest  release  dates  r,  j  of  the  operations,  then: 

/Rm  +  Ej  Ei  7 i,j  ' 


Crt  >  E 


M 


Proof:  This  demonstration  is  inspired  of  [9],  Tn  1987,  Carlier  have  proposed  an  interesting  lower  bound  for  the 
problem  of  n  operations  on  m  identical  machines  in  minimizing  the  makespan.  This  bound  takes  account  of  the  release 
date  of  each  operation  and  shows  that  the  activity  intervals  of  the  machines,  in  the  best  case,  have  the  following  form: 
[rk-  >  /*]  V  1  <  k  <  m  where  f*  is  the  optimal  value  of  the  makespan  and  [rq  ,  rj, , r, ;m  ]  the  m  smallest  release  dates  of 
the  operations.  By  using  the  same  idea  in  our  problem,  we  obtain  the  following  relation:  Crj  >  r 3  ,  Moreover,  we 

have  Cr 3  >  )Tb  "  ,  j  and  the  makespan  is  integer,  therefore,  the  lemma  is  justified.  ■ 

Definition  4:  Let  IN ,  the  set  of  the  combinations  const ituded  of  Nl  operations  among  the  NT  operations  and  Cjv, 
an  element  of  Eni-  CNi  =  {Oiujl,Oi2ij2,  We  define  7Lc;v,  as  follows: 


D-,cn,  =  ,  min 

1  <q<N 


1  <q<Nl 


Proposition  1 :  is  a  lower  bound  of  the  value  of  the  date  at  which  the  machine  Mk  can  carry  out  Nr  operations. 

Tk,Ni=  min  {'lk:cx, } 

Cn,£En, 

Proof:  For  a  given  combination  C^i  of  E^h  the  machine  Mk  can  start  to  perform  operations  at  least  at  the  date 
t  =  mini<9<jv/  {riq,jq  }  and  must  work  for  a  duration  at  least  equal  to  y^1<9< N,  diq,jq,k-  B 


4 


Lemma  S:  Let  Vz  the  subset  of  the  operations  set  defined  as  follows:  Vz  =  {Oi.+llj.+1 ,  Oj.+2j.+2ip,.,  OiNT,jNT  ]  such 
that  rit+ujt+1  <  ?’i.+2,j.+2  <  -  <  nNT,jNT,  then: 


Tk  Ni=  min  /',  .  I  <h 

'  l<z<NT-N>+ 1  \ 


+  ,  min 

C!  „y£ 


(a*  (c;iA„)) 


JV/  ^  a,iV/ 


where:  Ak  (C'z  Nl )  is  the  sum  of  the  processing  times  of  the  operations  of  C'  V/  on  M*;  C'  Y,  is  an  element  of  El  Nl 
and  El  Nl  is  the  set  of  the  combinations  of  (Nf  —  1)  operations  chosen  among  the  (NT  —  z)  operations  of  Vz  for 

Proof:  For  proof,  see  [10].  ■ 


Let  us  suppose  that  a  (k,  Nt)  =  1,  then  the  machine  Mk  will  work  until  a  date  at  least  equal  to  Tk<y,.  In  this  case, 
the  remainder  of  the  operations  (the  NT  —  Nl  operations)  must  be  carried  out  on  lJk .  This  obligation  can  be  regarded 
as  a  scheduling  problem  of  (NT  —  Nf)  operations  on  the  (M  —  1)  machines  of  Uk .  Therefore,  for  such  a  problem,  we 
can  apply  the  result  of  lemma  7  to  undervalue  the  necessary  time  for  such  an  execution  (the  minorant  of  such  necessary 

time  is  noted  A k,Ni)-  Thus,  if  NT  —  Nf  >  M  —  1,  A^.y,  =  E  m-i''  )  where  Rm- i  is  the  sum  of  the  (M  —  1) 

smallest  values  of  r, j.  Tn  the  contrary  case  and  if  NT  —  Nf  <  M  —  1 ,  it  is  more  interesting  to  undervalue  the  time  in 
question  by  E  ^  RNTNf'C\l!'N'  ^  where  Rnt  -Ni  is  the  sum  of  the  NT  —  Nf  smallest  values  of  r,  j .  Therfore,  we  have  the 
following  lemma. 

Lemma  9:  If  a  (k,  Nf)  =  1 ,  then  we  need  to  perform  operations  until  a  date  at  least  equal  to  T/,.  y,  such  that: 

Tk,N!  =  max  (Tk,Nh^k,Ni) 

with: 

A  Ay Vi  =  I  if  NT  ~  N'  ~  M  ~  1  - 

Aa .n,  =  e( RxT-x,  +  Tk  x'\  if  NT  _  N,  <  M  _  , 

V  NT-Nf  ) 

Theorem  1 :  rko  N  is  a  lower  bound  of  (Jr  | : 

Ti-n  n  =  min  \  min  {rk  Ni} 
ko’N  l<k<M  (Nl>N 

Proof:  min  V(>  y  {ta^y/}  represent  a  lower  bound  of  the  makespan  if  the  machine  Mk  carries  out  a  number  of 

operations  superior  or  equal  to  N.  Tn  addition,  there  is  at  least  one  machine  that  verify  such  a  condition,  therefore,  the 
theorem  is  justified.  ■ 

Let  1 7  the  set  of  all  the  operations  classified  in  the  ascending  order  according  to  the  values  of  rl:j  : 

I  —  {Oi1  ,ji  1  Oi-2,32  1  ■■■!  OiNTjNT  }  such  that  Ti1j1  <  I  i2j  2  <  ...  <  fiNTJNT 

Tt  is  clear  that  any  lower  bound  of  the  scheduling  problem  of  Vz  is  also  a  lower  bound  for  the  initial  problem  (the 
scheduling  of  V). 


Lemma  10:  LB2  is  an  improvement  of  the  lower  bound  proposed  in  lemma  7: 


LB  2  =  max 

I  <q<NT-M 


(  1  /h=q+M-i  h=mr  \  | 

M'  L  1  L  '  ( 

(  ^  \  h=q  h  =  q  )  ) 


Proof:  Such  a  result  is  justified  by  applying  lemma  7  to  the  subsets  I  Tz  for  1  <j<  NT  —  M  —  1 . 
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Corollary  1  :  Using  the  preceding  theorems  and  the  preceding  lemmas,  we  obtain  the  fo lowing  relation: 


Cri  >  max  (j/  ■  V  -  •  J  ,  LB2, Tko  ^j 


Remark 1  :  Tn  a  previous  work  [13],  we  have  showed  that  it  exists  an  equivalence  between  a  flexible  job-shop  with 
release  dates  and  a  flexible  job-shop  without  release  dates.  Thus,  we  can  use  this  equivalence  to  find  other  relations  as 
the  following  relation: 


Cr  i  >  max 


(J( 


Ej  tj  +  Ej  E,  "a,. 

N  +  M 


■  kn,Nl 


with  N/  =  E  ^  E+E  J '  E0  N i  =  mini<fc<M  \Dk  and  Dk  the  function  previously  defined  in  the  preceding  subsec¬ 
tion,  but  one  can  easily  show  that  such  bounds  are  less  interesting  than  those  defined  by  the  preceding  corollary. 


D.  Recapitulation 

The  preceding  minorations  will  enable  us  to  compute  some  limits  for  the  values  corresponding  to  the  three  criteria 
considered.  These  limits  are  defined  by  the  following  relations: 


Cri  >  Cr{ ,  Ct'2  >  Ci'2  and  Cr3  >  Cr| 

with: 


Cri  =  max 


(max 


,LB  2 


>  kn,N  ’ 


Ct’2  =  max 


Ej  E  /  tj 

M 


and  Cr3  =  SSv, 

j  * 


E.  Complexity 

The  formulas  of  the  different  lower  bounds  are  implemented  according  to  the  corresponding  algorithms.  All  these 
algorithms  are  polynomial.  Theirs  algorithmic  complexities  are  presented  in  Table  T. 

TV.  Simulation  Results 

To  test  the  efficiency  of  the  lower  bounds  in  the  evaluation  of  the  system  performance,  many  computational  experiments 
have  been  carried  out.  This  test  consists  in  applying  a  Controlled  Evolutionary  Approach  (CEA)  [2].  The  objective 
of  such  simulations  is  not  to  evaluate  the  efficiency  of  the  CEA  (in  fact,  the  performance  of  such  a  method  has  been 
demonstrated  in  previous  publications  [2]  [14]).  The  objective  considered  is  to  measure  the  quality  of  the  lower  bounds 
proposed  by  comparing  them  to  the  various  values  of  the  criteria  associated  to  the  solutions  given  by  the  evoked 
method.  Tn  this  section,  we  give  a  short  description  of  the  CEA,  then,  we  present  the  considered  examples  and  we  finish 
by  comparing  the  criteria  values  of  the  obtained  solutions  to  the  lower  bounds  values. 


TABLE  I 

Lower  bounds  complexities 


Lower  bound 

Complexity 

Crj 

O  (NT2) 

Cr:2 

O  (NT2) 

Cri 

O  (NT) 
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A.  Used  Method 

The  problem  considered  presents  two  main  difficulties.  The  first  one  is  the  assignment  of  each  operation  to  the 
suitable  machine.  The  second  difficulty  is  the  calculation  of  the  starting  times  t,  j  of  each  operation  O,  j.  To  solve  such 
a  problem,  we  apply  initially  an  Approach  by  Localization  [2] .  This  approach  is  a  heuristic  which  makes  it  possible 
to  assign  the  operations  to  the  machines  by  taking  account  of  the  processing  times  and  the  workloads  of  the  machines 
to  which  we  have  already  assigned  operations.  Then,  it  makes  it  possible  to  solve  the  problem  of  the  tasks  sequencing 
thanks  to  an  algorithm  called  “Scheduling  Algorithm”  [2]  which  calculates  the  starting  times  t\  j  by  taking  into  account 
the  availabilities  of  the  machines  and  the  precedence  constraints.  The  conflicts  are  solved  by  using  traditional  priority 
rules  (SPT,  LPT,  FTFO,  LTFO,  FTRO  [15],  [16]),  thus,  we  obtain  a  set  of  schedules  according  to  the  used  priority  rules. 
To  such  a  set,  we  apply  an  evolutionary  approach  which  is  based  on  the  schemata  theorem  introduced  in  the  genetic 
algorithms  field.  Such  an  approach  consists  in  the  design  of  an  assignment  model  which  will  be  useful  to  construct  the 
set  of  the  new  individuals.  The  objective  is  to  integrate  the  good  qualities  contained  in  the  schemata  [2],  [17],  [18]  in 
order  to  make  the  evolutionary  algorithm  more  effecient  and  more  rapid.  Tn  fact,  the  construction  of  the  solutions  is 
done  by  giving  the  priority  to  the  reproduction  of  the  individuals  respecting  the  model  generated  by  the  assignment 
schemata  and  not  starting  from  the  whole  set  of  the  chromosomes  (for  further  details,  the  reader  is  invited  to  consult 
[2]  and  [18]).  The  multiobjective  evaluation  of  the  solutions  is  carried  out  using  a  fuzzy  Pareto  approach  [14].  Such  an 
approach  is  based  on  the  weighted  aggregation  of  the  different  objective  functions  at  each  iteration  of  the  evolutionary 
algorithm.  The  particularity  of  such  an  approach  consists  in  the  fuzzy  computation  of  weights  by  giving  the  priority  to 
the  objective  functions  which  the  values  are  far  from  the  corresponding  lower  bound  value  [14]. 

R.  Results 

Many  series  of  examples  have  been  tested  to  evaluate  the  quality  of  the  lower  bounds  based  on  practical  data.  As 
an  example,  the  reader  could  consult  the  simulation  results  of  some  instances  at  the  web  address:  http://www.ec- 
lille.fr/~kacem/testsPareto.pdf.  These  instances  come  from  the  literature  [4],  [14]  and  present  problems  with  4  to  25 
jobs,  generally  using  10  machines  with  12  to  75  operations  with  total  flexibility.  The  various  results  are  summarized  in 
Table  IT.  For  each  instance,  we  present  the  values  of  the  different  lower  bounds  and  the  values  of  the  criteria  for  the 
Pareto  optimal  solutions  obtained  by  our  fuzzy  evolutionary  approach. 

The  results  presented  in  Table  TT  show  that  the  solutions  obtained  are  generally  very  close  to  the  optimal  one.  Tn  the 
case  of  the  Parallel  Machines  problem  (a  particular  case  of  the  flexible  job-shop  described  by  the  instances  7g,  /g,  I7.  Tg , 
7g  and  7io),  we  obtain  no  distance  between  the  lower  bounds  values  and  the  different  obtained  solutions.  Tn  the  general 
case,  the  small  distance  that  we  can  have  is  due  to  the  difficulty  of  multiobjective  optimization  which  considers  several 
nonhomogeneous  and  antagonistic  criteria  at  the  same  time  [19].  Recently,  we  have  also  tested  some  instances  coming 
from  the  benchmarks  of  Hurink  [20].  The  results  obtained  confirm  the  good  quality  of  the  proposed  lower  bounds. 

V.  Conclusion 

Tn  this  paper,  \\7e  have  proposed  a  set  of  lower  bounds  to  make  it  easier  the  multiobjective  evaluation  problem  of  a 
schedule  performance  in  the  case  of  flexible  job  shops.  These  lower  bounds  make  it  possible  to  estimate  precise  limits 
for  the  optimal  values  of  the  corresponding  criteria.  Tn  fact,  different  simulations  show  that  the  little  distance  between 
such  limits  and  the  values  of  the  criteria  obtained  for  the  solutions  generated  by  the  evolutionary  fuzzy  approach  is 
generally  satisfactory  and  promising.  This  result  is  very  interesting  to  facilitate  the  study  of  others  multiobjective 
concepts  (Uniform  Design  concept  [19])  that  we  will  consider  as  perspective  in  our  future  work. 
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Instances 
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